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Let / be a compact d-dimensional manifold, let X: 7 —> 7?. be a 
Gaussian process with regular paths and let Fi{u), u €TZ, be the 
probability distribution function of supjgjA'(t). 

We prove that under certain regularity and nondegeneracy condi¬ 
tions, Fi is a C^-function and satisfies a certain implicit equation that 
permits to give bounds for its values and to compute its asymptotic 
behavior as u—> -|-oo. This is a partial extension of previous results 
by the authors in the case d=l. 

Our methods use strongly the so-called Rice formulae for the mo¬ 
ments of the number of roots of an equation of the form Z{t) = x, 
where Z\I^ TiS is a random field and a: is a hxed point in TiS. 

We also give proofs for this kind of formulae, which have their own 
interest beyond the present application. 

1. Introduction and notation. Let / be a d-dimensional compact man¬ 
ifold and let A: I —> 7^ be a Gaussian process with regular paths defined 
on some probability space (n,A.,P). Define Mj = siip^^jX{t) and F[(u) = 
P{M/ <u}, u€TZ, the probability distribution function of the random vari¬ 
able Mj. Our aim is to study the regularity of the function Fj when d> 1. 

There exist a certain number of general results on this subject, starting 
from the papers by Ylvisaker (1968) and Tsirelson (1975) [see also Weber 
(1985), Lifshits (1995), Diebolt and Posse (1996) and references therein]. The 
main purpose of this paper is to extend to d > 1 some of the results about the 
regularity of the function u Fj{u) in Azai's and Wschebor (2001), which 
concern the case d = l. 

Our main tool here is the Rice formula for the moments of the number of 
roots N^{I) of the equation Z(t) =u on the set I, where {Z(t): t G /} is an 
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7^'^-valued Gaussian field, / is a subset of and u is a given point in TZ'^. 
For d > 1, even though it has been used in various contexts, as far as the 
authors know, a full proof of the Rice formula for the moments of N^{I) 
seems to have only been published by Adler (1981) for the first moment of 
the number of critical points of a real-valued stationary Gaussian process 
with a d-dimensional parameter, and extended by Azai's and Delmas (2002) 
to the case of processes with constant variance. Gabaha (1985) contains 
related formulae for random fields; see also the Ph.D. thesis of Konakov 
cited by Piterbarg (1996b). In the next section we give a more general result 
which has an interest that goes beyond the application of the present paper. 
At the same time the proof appears to be simpler than previous ones. We 
have also included the proof of the formula for higher moments, which in 
fact follows easily from the first moment. Both extend with no difficulties to 
certain classes of non-Gaussian processes. 

It should be pointed out that the validity of the Rice formula for Lebesgue- 
almost every u G TZ'^ is easy to prove [Brillinger (1972)] but this is insufficient 
for a certain number of standard applications. For example, assume A : / ^ 
A is a real-valued random process and one is willing to compute the moments 
of the number of critical points of X. Then, we must take for Z the random 
field Z[t) = X'{t) and the formula one needs is for the precise value u = 0 
so that a formula for almost every u does not solve the problem. 

We have added the Rice formula for processes defined on smooth mani¬ 
folds. Even though the Rice formula is local, this is convenient for various 
applications. We will need a formula of this sort to state and prove the im¬ 
plicit formulae for the derivatives of the distribution of the maximum (see 
Section 3). 

The results on the differentiation of Fj are partial extensions of Azais and 
Wschebor (2001). Here, we have only considered the first derivative Fj(u). 
In fact, one can push our procedure one step more and prove the existence 
of Fj{u) as well as some implicit formula for it. But we have not included 
this in the present paper since formulae become very complicated and it 
is unclear at present whether the actual computations can be performed, 
even in simple examples. The technical reason for this is that, following the 
present method, to compute F"{u), one needs to differentiate expressions 
that contain the “helix process” that we introduce in Section 4, containing 
singularities with unpleasant behavior [see Azais and Wschebor (2002)]. 

For Gaussian fields defined on a d-dimensional regular manifold (d > 1) 
and possessing regular paths we obtain some improvements with respect to 
classical and general results due to Tsirelson (1975) for Gaussian sequences. 
An example is Gorollary 5.1, which provides an asymptotic formula for Fj{u) 
as tt —> -|-oo which is explicit in terms of the covariance of the process and 
can be compared with Theorem 4 in Tsirelson (1975) where an implicit 
expression depending on the function F itself is given. 
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We use the following notation: 

If Z is a smooth function U TZ'^ , U a subset of TZ^, its successive 
derivatives are denoted Z', Z" ,..., and considered, respectively, as lin¬ 
ear, bilinear, ..., fe-linear forms on 7Z^. For example, X^^\t)[vi,V 2 ,V 3 ] = 

J2i,j,k=i . The same notation is used for a derivative on a C°° 

manifold. 

/, dl and I are, respectively, the interior, the boundary and the closure 
of the set I. If ^ is a random vector with values in TZ'^, whenever they exist, 
we denote by p^{x) the value of the density of ^ at the point x, by E(^) 
its expectation and by Var(^) its variance-covariance matrix. A is Lebesgue 
measure. If u,v are points in TZ^, {u,v) denotes their usual scalar product 
and ||tt|| the Euclidean norm of u. Eor M a d x d real matrix, we denote 
ll^-^ll =sup||^ll=i ||Mx||. 

Also for symmetric M, M y 0 (resp. M ~<0) denotes that M is positive 
dehnite (resp. negative definite). denotes the complement of the set A. 
Eor real x, = sup(x,0), x~ = sup(—x,0). 

2. Rice formulae. Our main results in this section are the following: 

Theorem 2.1. Let Z:I I a compact subset oflZ^, be a random 

field and u G TZ'^. Assume that: 

(AO) Z is Gaussian. 

(Al) t'^Z{t) is a.s. of class . 

(A2) For each t £ I, Z{t) has a nondegenerate distribution [i.e., Var(Z(t)) 

0]- 

(A3) F{3t G I, Z{t) = u, det{Z'{t)) = 0} = 0. 

(A4) A(dl) = 0. 

Then 

(1) E(iVf(/)) = jF{\det{Z'm/Z{t)=u)pzi^t){u)dL 

and both members are finite. 

Theorem 2.2. Let k, k>2, be an integer. Assume the same hypotheses 
as in Theorem 2.1 except for (A2), which is replaced by: 

(A'2) for ti,... ,tk £ I pairwise different values of the parameter, the distri¬ 
bution of {Z{ti),..., Z{tk)) does not degenerate in Then 

F[{N^{L)){N^{I) - 1) • • • - A: + 1)] 

( k 

n I det(Z'(t,))|/Z(fi) = • • • = Z(4) = u 

3 = 1 

X Pz{ti),...,z{tt) (..u,...,u)dti--- dtk, 
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where both members may be infinite. 

Remark. Note that Theorem 2.1 (resp. Theorem 2.2) remains valid if 
one replaces I by / in (1) or (2) and if hypotheses (A0)-(A2) [resp. (A'2)] 
and (A3) are verified. This follows immediately from the above statements. 
A standard extension argument shows that (1) holds true if one replaces I 
by any Borel subset of I. 

Sufficient conditions for hypothesis (A3) to hold are given by the next 
proposition. Under condition (a) the result is proved in Lemma 5 of Cucker 
and Wschebor (2003). Under condition (b) the proof is straightforward. 

Proposition 2.1. Let Z:I I a compact subset ofTZ'^, be a ran¬ 

dom field with paths of class and u G TZ'^. Assume that: 

(i) Pz{t){x) < C for all t £ I and x in some neighborhood of u. 

(ii) At least one of the two following hypotheses is satisfied: 

(a) a.s. t ^ Z{t) is of class , 

(b) a{5) = supjg^ P{| det(Z'(t))| < 6/Z{t) = x} ^ 0 as <5 ^ 0, 

where V{u) is some neighborhood of u. 

Then (A3) holds true. 

The following lemma is easy to prove. 

Lemma 2.1. With the notation of Theorem 2.1, suppose that (Al) and 
(A4) hold true and that Pz(t) (^) ^ C for allt G I and x in some neighborhood 
ofu. Then V{N^{dI) / 0} = 0. 

Lemma 2.2. Let Z '.L ^ TZ^, L a compact subset ofTZ'^, be a function 
and u a point in TZ'^. Assume that: 

(a) inftg 2 -i(W)(^min(^'(i))) > A > 0, 

(b) uz'ig) < A/d, 

where loz' is the continuity modulus of Z', defined as the maximum of the 
continuity moduli of its entries, Amin(Af) is the square root of the smallest 
eigenvalue of M'^M and rj is a positive number. 

Then, ifti,t 2 are two distinct roots of the equation Z{t) =u such that the 
segment [ti,t 2 ] is contained in L, the Euclidean distance between ti and t 2 
is greater than rj. 
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Proof. Set fj = ||ti — t 2 ||, v = | | *^_*^|| . Using the mean value theorem, 
for z = 1,..., d, there exists G [^ 1 ,^ 2 ] such that {Z'(^i)v)i = 0. Thus 

d d 

< \Z\ti)ik - Z\(i)ik\\vk\ <uJz'{fi) Y l^fcl - ^Z'{fi)^fd- 

k=l k=l 

In conclusion, A < \^\a{Z' {ti)) < ||Z'(ti)u|| < iOz'{fj)d, which implies fjy rj. 

□ 

Proof of Theorem 2.1. Consider a continuous nondecreasing func¬ 
tion F such that F{x) = 0 for x < 1/2, F{x) = 1 for a; > 1. Let A and rj be 
positive real numbers. Define the random function 

aA,riiu)=F(^-^mi[X^iniZ'{s)) + \\Z{s)-u\\]^ x (^1-f(^^u;z'{v)^'^, 

(3) 

and the set = {t € I :\\t — s|| > ry,Vs ^ /}. If aA,ri{u) > 0 and N^{I^ri) 
does not vanish, conditions (a) and (b) in Lemma 2.2 are satished. Hence, 
in each ball with diameter ^ centered at a point in there is at most 
one root of the equation Z{t) = u, and a compactness argument shows that 
is bounded by a constant I), depending only on r] and on the 

set /. 

Take now any real-valued nonrandom continuous function /: —> TZ 
with compact support. Because of the coarea formula [Federer (1969), The¬ 
orem 3.2.3], since a.s. Z is Lipschitz and aA,r){y) ■ f{u) is integrable, 

/ f{u)N^{I-r,)aA,r,{u)du= [ \det{Z'{t))\f{Z{t))aA,rj{Z{t))dt. 

Jl-r, 

Taking expectations in both sides, 

/ f{u)E{N^{I_rj)aA,rj{u))du 

J-Rd- 

= [ fiu)du f E{\det{Z'{t))\aA,r,iu)/Z{t) = u)pz(t)iu)dt. 

J J 1—r) 

It follows that the two functions 

(i) E{N^{I_^)aA,r,{u)), 

(ii) !i_^^{\dQ^{Z'{t))\aA,r,{u)/Z{t) = u)pz{t){u)dt, 

coincide Lebesgue-almost everywhere as functions of u. 

Let us prove that both functions are continuous, hence they are equal for 
every u G TZ'^. 
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Fix It = Wo and let us show that the function in (i) is continuous at w = wq- 
Consider the random variable inside the expectation sign in (i). Almost 
surely, there is no point t in Z~^{{uo}) such that det{Z'{t)) = 0. By the lo¬ 
cal inversion theorem, Z{-) is invertible in some neighborhood of each point 
belonging to Z~^{{uo}) and the distance from Z(t) to uq is bounded below 
by a positive number for t G I-rj outside of the union of these neighborhoods. 
This implies that, a.s., as a function of u, is constant in some (ran¬ 

dom) neighborhood of uq. On the other hand, it is clear from its definition 
that the function u ^ is continuous and bounded. We may now ap¬ 

ply dominated convergence as w —> wq, since {I-rj)aA,r]{'U‘) is bounded by 
a constant that does not depend on u. 

For the continuity of (ii), it is enough to prove that, for each t € I the 
conditional expectation in the integrand is a continuous function of u. Note 
that the random variable \det{Z'{t))\a^^riiu) is a functional defined on 
{{Z(s), Z'(s)): s £ 1}. Perform a Gaussian regression of {Z{s), Z'{s)): s G I 
with respect to the random variable Z{t), that is, write 

Z{s) = Y\s) + a\s)Z{t), 

Z'{s) = Ylis)+f3j{s)Z{t), j = 

where Zj{s), j = 1,... ,d, denote the columns of Z'{s), T*(s) and Yj{s) are 
Gaussian vectors, independent of Z{t) for each s G /, and the regression 
matrices Q:*(s), Pj{s), j = 1,... ,d, are continuous functions of s,t [take into 
account (A2)]. Replacing in the conditional expectation, we are now able 
to get rid of the conditioning, and using the fact that the moments of the 
supremum of an a.s. bounded Gaussian process are finite, the continuity in 
u follows by dominated convergence. 

So, now we fix w G and make ?7 j 0, A | 0 in that order, both in 
(i) and (ii). For (i) one can use Beppo Levi’s theorem. Note that almost 
surely [I-r]) T dY^{I) = (I), where the last equality follows from Lemma 

2.1. On the other hand, the same Lemma 2.1 plus (A3) imply together that, 
almost surely, 

inf[Amm(^'(s)) + ll^('S) - m||] > 0 

S&I 

so that the first factor in the right-hand side of (3) increases to 1 as A 
decreases to zero. Hence by Beppo Levi’s theorem, 

Jim E ((/_^) a A,r, (w)) = E( (I)). 

For (ii), one can proceed in a similar way after deconditioning obtaining (1). 
To finish the proof, remark that standard Gaussian calculations show the 
finiteness of the right-hand side of (1). □ 
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Proof of Theorem 2.2. For each 5 > 0, define the domain 
Dk,5{I) = {(ii, - ■ • ,4) e \\ti - tjW > (5 if i / j,z, j = 1,... ,/c} 
and the process Z 

(ti,... ,4) S Dk^sil) Z{ti,... ,tk) = {Z(ti). ,Z{tk)). 

It is clear that Z satishes the hypotheses of Theorem 2.1 for every value 
So, 

(4) = / E(nidet(Z'(4))|/Z(4) = ... = Z(4) = u') 

Jd^Ai) \f=i ) 

X Pz{ti),...,z(tk){u, ■ ■ ■ ,u) dti ■ ■ ■ dtk. 

To hnish, let 5 i 0, note that {N^{!) — 1)... {N^{I) — A: + 1) is 
the monotone limit of and that the diagonal Dk{I) = 

{(ti,..., 4) ^ ti = tj for some pair i,j, i / j} has zero Lebesgue measure 
in {TZ^f. 

□ 

Remark. Even thought we will not use this in the present paper, we 
point out that it is easy to adapt the proofs of Theorems 2.1 and 2.2 to 
certain classes of non-Gaussian processes. 

Eor example, the statement of Theorem 2.1 remains valid if one replaces 
hypotheses (AO) and (A2), respectively, by the following (BO) and (B2): 

(BO) Z{t) = H{Y{t)) for t€ I, where Y: I ^ is a Gaussian process with 
paths such that for each t G I, Y(t) has a nondegenerate distribution 
and H: is a function. 

(B2) For each t G I, Z{t) has a density Pz{t) which is continuous as a function 
of {t,u). 

Note that (BO) and (B2) together imply that n>d. The only change to 
be introduced in the proof of the theorem is in the continuity of (ii) where 
the regression is performed on Y(t) instead of Z{t). 

Similarly, the statement of Theorem 2.2 remains valid if we replace (AO) 
by (BO) and add the requirement that the joint density of Z{ti ),..., Z(tk) 
be a continuous function of ti,..., 4; 'U for pairwise different ti,... ,tk- 
Now consider a process X from I to TZ and define 

= (j{t G I,X{-) has a local maximum at the point t,X{t) > u}, 
M^^I) = e I, X'{t) = 0, X{t) > u}. 
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The problem of writing Rice formulae for the factorial moments of these 
random variables can be considered as a particular case of the previous one 
and the proofs are the same, mutatis mutandis. For further use, we state 
as a theorem the Rice formula for the expectation. For breavity we do not 
state the equivalent of Theorem 2.2, which holds true similarly. 

Theorem 2.3. Let X ■. I TZ, I a compact subset ofTZ'^, be a random 

field. Let u^TZ, define i = 1,2, as above. For each dx d real sym¬ 
metric matrix M, we put 5^{M) := | det(M)|]lM-< 0 ) := |det(M)|. 

Assume: 

(AO) X is Gaussian, 

(A"l) a.s. t'^X{t) is of class , 

{X"2) for each t £ I, X(t),X'{t) has a nondegenerate distribution in TZ^ x 

n<^, 

(A"3) either a.s. tX{t) is of class or a((5) = supjgj 3 ,/gy(Q) P(| det(A"(t))| 
5/X'{t) = x') ^ 0 as 6 ^ 0, where F(0) denotes some neighborhood 
ofO, 

(A4) dl has zero Lebesgue measure. 

Then, for i = l,2, 

poo p 

E(M„^,(/)) = / dx E{dfiX''it))/Xit) =x,X'{t) = 0)px(t),x'it){x,0)dt 
Ju JI 

and both members are finite. 

2.1. Processes defined on a smooth manifold. Let 1/ be a differentiable 
manifold (by differentiable we mean infinitely differentiable) of dimension d. 
We suppose that U is orientable in the sense that there exists a nonvanishing 
differentiable d-form LI on U. This is equivalent to assuming that there exists 
an atlas {{Ui,f>i);i G I) such that for any pair of intersecting charts 
{Lfj^fij), the Jacobian of the map fii o is positive. 

We consider a Gaussian stochastic process with real values and paths 
X = {A(t) :t gU} defined on the manifold U. In this section we first write 
Rice formulae for this kind of processes without further hypotheses on U. 
When U is equipped with a Riemannian metric, we give, without details and 
proof, a nicer form. Other forms exist also when U is naturally embedded 
in a Euclidean space, but we do not need this in the sequel [see Azai's and 
Wschebor (2002)]. 

We will assume that in every chart X{t) and DX{t) have a nondegenerate 
joint distribution and that hypothesis (A''3) is verified. Eor S a Borel subset 
of U, the following quantities are well defined and measurable: the 

number of local maxima and M^ 2 {S), the number of critical points. 
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Proposition 2.2. For k = l,2 the quantity which is expressed in every 
chart (j) with coordinates si,..., as 


(5) 


+ 00 


dxF{5^{Y''{s))/Y{s) = x,Y'{s) = t))pY(s),Y'{s){x,Q) f\ dsi, 


i=l 

k-l 


where l^(s) is the process X written in the chart: Y = X o cj) defines a 
d-form on U and for every Borel set S <ZU, 




Proof. Note that a d-form is a measure on U whose image in each 
chart is absolutely continuous with respect to Lebesgue measure f\f=idsi. 
To prove that (5) defines a d-form, it is sufficient to prove that its density 
with respect to Af=i dsi satisfies locally the change-of-variable formula. Let 
{U 2 , 4 > 2 ) be two intersecting charts and set 

U3-.= UinU 2 ; Yi:=Xocj)f^; Y 2 := X o ■ H := (j )2 o 
Denote by sj and sf, i = i,..., d, the coordinates in each chart. We have 


dYi 

dsi 

d^Yi 
dsj dsj 


E 

E 


dY2 dHij 


dsi 

d‘^Y2 dHi, dHj, 
ds] ds] 


^ dY2 d^H,, 
^ dsi ' 


Thus at every point 

Yl{s^) = {H'{s^)fYf{s^), 


Pyi(M),y/(M)(a:,0) = pY2{s2),Y'is^){x,0)\ det{H'{s^))\ \ 
and at a singular point, 

Yfi{s^) = {H'{s^)fYfi{s^)H\s^). 

On the other hand, by the change-of-variable formula, 
d d 

/\ ds,^ = I det(id'('S^))r^ A ^'5?' 

i=l i=l 


Replacing in the integrand in (5), one checks the desired result. 

For the second part again it suffices to prove it locally for an open subset 
S included in a unique chart. Let (S', fi) be a chart and let again Y (s) be 
the process written in this chart. It suffices to check that 


E(M,A(S)) 

(b) I" r +00 

= dX{s) dxE{5^{Y"{s))/Y{s) =x,Y'{s) = 0)pY(s),Y'is)ix,0)- 

J 4>{S) J u 
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Since is equal to M^^{(/)(5)}, we see that the result is a direct 

consequence of Theorem 2.3. 

Even though in the integrand in (5) the product does not depend on the 
parameterization, each factor does. When the manifold U is equipped with 
a Riemannian metric it is possible to rewrite (5) as 

r+oo 

(7) / dx^{5\V^X{s))/X{s)=x,VX{s) = {))pxis),vx{s){x,^)\o\, 

Ju 

where V^X(s) and VAl(s) are respectively the Hessian and the gradient 
read in an orthonormal basis. This formula is close to a formula by Taylor 
and Adler (2002) for the expected Euler characteristic. 

Remark. One can consider a number of variants of Rice formulae, in 
which we may be interested in computing the moments of the number of 
roots of the equation Z{t) =u under some additional conditions. This has 
been the case in the statement of Theorem 2.3 in which we have given 
formulae for the first moment of the number of zeroes of X' in which X is 
bigger than u (i = 2) and also the real-valued process X has a local maximum 
(i = l). 

We just consider below two additional examples of variants that we state 
here for further reference. We limit the statements to random fields dehned 
on subsets of 77'^. Similar statements hold true when the parameter set is 
a general smooth manifold. Proofs are essentially the same as the previous 
ones. 

Variant 1. Assume that Z\ , Z 2 are TT.'^-valued random fields defined on 
compact subsets /i ,/2 of 77'^ and suppose that {Zi,Ii), f = 1,2, satisfy the 
hypotheses of Theorem 2.1 and that for every s G R and t G I 2 , the distribu¬ 
tion of (Zi(s), ^ 2 ( 7 )) does not degenerate. Then, for each pair ui,U 2 G 77'^, 

E(<nAX^(^2)) 

( 8 ) =[ dtidt2E{\det{Z[{ti))\\det{Z2{t2))\/Zi{ti) =ui,Z2it2) =U2) 

Jhxh 

X PZi(q),Z2(t2)(^D“2)- 

Variant 2. Let Z,I be as in Theorem 2.1 and let be a real-valued 
bounded random variable which is measurable with respect to the cr-algebra 
generated by the process Z. Assume that for each 7 G I, there exists a contin¬ 
uous Gaussian process {V*(s): s G /}, for each s, 7 G I a nonrandom function 
a*(s): 77'^ ^ IZ'^ and a Borel-measurable function g:C ^IZ where C is space 
of real-valued continuous functions on I equipped with the supremum norm, 
such that: 


DISTRIBUTION OF THE MAXIMUM 


11 


1. C = 9{Y\-)+a\-)Z{t)), 

2. Y^{-) and Z{t) are independent, 

3. for each uq G TZ, almost surely the function u ^ g{Y^{-) + a^{-)u) is con¬ 
tinuous at u = uo. 

Then the formula 

E{N^{m = jE{\detiZ\tm/Z{t)=u)pzit){u)dt 

holds true. 

We will be particularly interested in the function ^ = 1mi<v for some 
V €TZ. We will see later on that it satisfies the above conditions under 
certain hypotheses on the process Z. 


3. First derivative, first form. Our main goals in this and the next sec¬ 
tion are to prove existence and regularity of the derivative of the function 
u ^ Fi{u) and, at the same time, that it satisfies some implicit formulae 
that can be used to provide bounds on it. In the following we assume that 
/ is a d-dimensional manifold embedded in TZ^^ N > d. a and a are 
respectively the geometric measures on I and dl. Unless explicit statement 
of the contrary is made, the topology on I will be the relative topology. 

In this section we prove formula (10) for Fj{u )—which we call “first 
form”—which is valid for A-almost every tt, under strong regularity con¬ 
ditions on the paths of the process X. In fact, the hypothesis that X is 
Gaussian is only used in the Rice formula itself and in Lemma 3.1 which 
gives a bound for the joint density Px{s),x{t),x'{s),x'{t)- fo both places, one 
can substitute Gaussianity by appropriate conditions that permit to obtain 
similar results. 

More generally, it is easy to see that inequality (9) is valid under quite 
general non-Gaussian conditions and implies that we have the following up¬ 
per bound for the density of the distribution of the random variable Mj: 


(9) 


F'j{u) < ^ E{5\X'^it))/Xit) = u, X\t) = 0)pxit),x'it) (u, O)a(dt) 

+ E(di(X"(t))/X(t) = n, W(t) = (u,O)d(dt), 


where the function has been dehned in the statement of Theorem 2.3 and 
X denotes the restriction of X to the boundary dF 

Even for d = l (one-parameter processes) and X Gaussian and stationary, 
inequality (9) provides reasonably good upper bounds for Fj{u) [see Diebolt 
and Posse (1996) and Azai's and Wschebor (2001)]. We will see an example 
for d = 2 at the end of this section. 

In the next section, we are able to prove that Ff{u) is a function 
and that formula (10) can be essentially simplified by getting rid of the 
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conditional expectation, thus obtaining the “second form” for the derivative. 
This is done under weaker regularity conditions but the assumption that X 
is Gaussian becomes essential. 


Definition 3.1. Let X: / —> 7?. be a real-valued stochastic process de¬ 
fined on a subset of We will say that X satisfies condition (Hk), k a 
positive integer, if the following three conditions hold true: 

(i) X is Gaussian; 

(ii) a.s. the paths of X are of class C^; 

(iii) for any choice of pairwise different values of the parameter ti,... ,tn, 

the joint distribution of the random variables X{ti),... ,X{tn), X'(ti),..., X'{tn ),..., X^^'^ (ti),..., X^^'^ ( 
has maximum rank. 


The next proposition shows that there exist processes that satisfy {Hk). 

Proposition 3.1. Let X = {X{t) -.t ^ TZ'^ he a centered stationary Gaus¬ 
sian process having continuous spectral density . Assume that f^{x) > 0 
for every x £ and that for any a > 0 f^{x) < CallxH"" holds true 

for some constant Ca and all x £ . Then, X satisfies {Hk) for every 
k = l,2,.... 


Proof. The proof is an adaptation of the proof of a related result for 
d=l [Cramer and Leadbetter (1967), page 203]; see Azai's and Wschebor 
( 2002 ). □ 


Theorem 3.1 (First derivative, first form). Let X :I be a Gaussian 
process, L a C°° compact d-dimensional manifold. Assume that X verifies 
(Hk) for every A: = 1,2,.... 

Then, the function u ^ Fi{u) is absolutely continuous and its Radon- 
Nikodym derivative is given for almost every u by 

F'j{u) = {-l)'^ jF{Aei{X"{t))lM,<u/X{t) = u,X'{t) = 0) 


( 10 ) 


^ PX{t),X'{t){u,^)cF{dt) 

+ (-1)'^-! [ E{dei{X"{t))lM,<u/X{t) = u, X'{t) = 0) 

Jdi 




Proof. For u <v and S (resp. S) a subset of I (resp. dl), let us denote 
Mu,v{S) =W{t£S:u< X{t) < V, X'{t) = 0, X''{t) -< 0}, 

Mu,v{S) = \){t£S:u<X{t)<v,X'{t) = 0,X''{t)-i0}. 
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Step 1. Let h> 0 and consider the increment 
Fj{u) — Fi{u — h) 

= P({M, < n} n [{M^-h,u{i) > 1} U {M^-h,um > I}])- 
Let us prove that 

(11) > 1,M„_,,^„(5/) > 1) = o(/i) as/ijO. 

In fact, for (5 > 0, 

> lA-h,u{dI) > 1 ) 

( 12 ) 


< F{Mu-h,u{.I-5)Mu-h,um) + E{Mu-h,u{I \ Is))- 

The first term in the right-hand side of (12) can be computed by means 
of a Rice-type formula, and it can be expressed as 



xdl 


a{dt)a{di) 


dx dx 


u—h 


X F[5^{X"{t))5^{X"{i))/X{t) = x,X{i) 


x,X'{t) 


0,X'(t) = 0) 


where the function 5^ has been defined in Theorem 2.3. 

Since in this integral ||t —1|| >S, the integrand is bounded and the integral 
is 0{h^). 

For the second term in (12) we apply the Rice formula again. Taking into 
account that the boundary of I is smooth and compact, we get 


E(M„_;,,„(I\/_5)) 

= / cr{dt) f E{5^{X''{t))/X{t) =x,X'{t) =0)px(t),x'it)ix,0)dx 

Jl\I-S Ju—h 

< (const)/i(T(/ \ I-s) < (const)M, 

where the constant does not depend on h and 5. Since 5 > 0 can be chosen 
arbitrarily small, (11) follows and we may write as /i —> 0: 

Fi (u) — Fj(u — h) 

= F{Mi < u, Mu-h,uii) > 1) + Mu-h,nidi) > 1) + oih). 

Note that the foregoing argument also implies that Fj is absolutely con¬ 
tinuous with respect to Lebesgue measure and that the density is bounded 
above by the right-hand side of (10). In fact, 

Fj{u) - Fi{u -h)< V{Mu-h,u{i) > 1) + ^{Mu-h,u{dl) > 1) 
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and it is enough to apply the Rice formula to each one of the expectations 
on the right-hand side. 

The delicate part of the proof consists in showing that we have equality 
in (10). 


Step 2. For g:I —>Tl we put H^Hoo = sup^^j |g'(t)| and if fe is a non¬ 
negative integer, \\g\\oc,k = supkT_+k2+-+ka<k\\dki,k2,...,ka9\\oo- For fixed 7 > 
0 (to be chosen later on) and h > 0,we denote by Eh = {||^||oo ,4 < h~'^}. 
Because of the Landau-Shepp-Fernique inequality [see Landau and Shepp 
(1970) or Fernique (1975)] there exist positive constants Ci,C 2 such that 

P{Eh) <Ciex.p[—C2h~'^'^]=o{h) as/i^O, 

so that to have (10) it suffices to show that, as /i —> 0, 

(13) 'E{[Mu-h,u{i) — = o{h), 

(14) E{[Mu-h,u{dI) - JlM/<ul£;h) = o{h). 

We prove (13). Equation (14) can be proved in a similar way. 

Put Mu-h,u = ^u-h,u{I)- We have, on applying the Rice formula for the 
second factorial moment, 


(15) 


E{\Mu-h,u - 


< E[Mu-h,u{Mu-h,u - ^)'^Eh) = As^ta{ds)a{dt), 

J Jlxl 


where 
As.t — 


dxi dx2 


u—h 


(16) 


X E(| det (X"(s)) det {X"{t))\lx/gs)^o,x"(t)^o'^Eh/^is) — 

Xit)=X 2 ,X'{s) = 0,X\t) = 0) 


X PX{s),X{t),X'{s),X'{t) {xiX2,0, 0). 

Our goal is to prove that Ag^t is o{h) as /i | 0 uniformly on s,t. Note that 
when s, t vary in a domain of the form Ds := {t,s € I: ||t — s|| > <5} for some 
5 > 0, then the Gaussian distribution in (16) is nondegenerate and Ag^t 
is bounded by (const)/i^, the constant depending on the minimum of the 
determinant: det Var(X(s), AL(t),X'(s),X'(t)), for s,t£Ds. 

So it is enough to prove that Ag^t = o{h) for ||t — s|| small, and we may 
assume that s and t are in the same chart {U, 4>)- Writing the process in this 
chart, we may assume that / is a ball or a half ball in 77'^. Let s, t be two such 
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points, and define the process V = F®’* by y(r) = X(s + T(t — s)); r e [0,1]. 
Under the conditioning one has 

y(o) = xi, r(i) = x2, Y'(o) = Y'(i) = o, 

Y"(0) = X"(s)l(t - s), (t - s)], Y"(l) = X"mt - s), (t - s)]. 

Consider the interpolation polynomial Q of degree 3 such that 
Q(0)=a:i, Q(1)=X2, Q'(0) = Q'(1) = 0. 

Check that 

Q(y) = xi + (x 2 - xi)y^(3 - 2y), (5"(0) = -Q"(l) = 6 (x 2 - xi). 

Denote Z{t) = Y{t) — Q(r),0 < r < 1. Under the conditioning, one has 
Z(0) = Z{1) = Z'{0) = Z'{\) = 0 and if also the event Eh occurs, an ele¬ 
mentary calculation shows that for 0 < r < 1 , 

(17) \Z''{t)\< sup sup (const)lit-sf/i-T-. 

re[o,i] 2 ! 2 ! 

On the other hand, check that if ^ is a positive semidehnite symmetric 
dx d real matrix and vi is a vector of Euclidean norm equal to 1, then the 
inequality 

(18) det( 74 ) < (^ui, ui) det(i?) 

holds true, where B is the {d—l)x{d— 1) matrix B = {{{Avj,Vk)))j^k= 2 ,...,d 
and {ui,U 2 ,... ,Vd} is an orthonormal basis of R'^ containing vi. 

Assume A"(s) is negative definite, and that the event Eh occurs. We 
can apply (18) to the matrix A = —X"{s) and the unit vector vi = {t — 
s)/||t — s||. Note that in that case, the elements of matrix B are of the form 
{—X"{s)vj,Vk), hence bounded by (const)h“'>'. So, 

det[-A"(s)] < {-X"{s)vi,vi)Cd = Cd[Y"{0)]-\\t-s\\-^h-^‘^-^^\ 

the constant depending only on the dimension d. 

Similarly, if X''{t) is negative definite, and the event Eh occurs, then 

det[-A"(t)] < Cd[Y''{l)]-\\t - s\\-^hAd-^h^ 

Hence, if C is the condition {A(s) = xi,X{t) = X 2 , X'{s) = 0, X'(t) = 0}, 

E(| det(A"(s))det(A"(t))|]ljf//(^)_.o,x"(t)-<olu^/<^) 

< Cj - s\\-^E{[Y"{0)]-[Y'' {1)]-1eJC) 

<Cj /i-2('='-i)T'||t-s||-4E(^ 

= Cj /l-2('^-l)7||i_s||-4E|^ 



< (const)Cj h — s 


4 
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We now turn to the density in (15) using the following lemma which is 
similar to Lemma 4.3, page 76, in Piterbarg (1996a). The proof is omitted. 

Lemma 3.1. ForaUs,t^I, 

(19) \\t - Px(s),x(L,N'(s),x'(t)(0,0,0,0) < D, 

where D is a constant. 

Back to the proof of the theorem, to bound the expression in (15) we use 
Lemma 3.1 and the bound on the conditional expectation, thus obtaining 

^{Mu-h,u{Mu-h,u - 1)1E;i) 

(20) < (const)C'J f f p — dsdt f f dxidx 2 

J Jlxl JJu—h 

< (const)/i2-2'='T' 

since the function {s,t) ^ p — is Lebesgue-integrable in / x /. The 

last constant depends only on the dimension d and the set I. Taking 7 small 
enough, (13) follows. □ 

Example. Let {AL(s,t)} be a real-valued two-parameter Gaussian, cen¬ 
tered stationary isotropic process with covariance T. Assume that r(0) = 

1 and that the spectral measure fj, is absolutely continuous with density 
^{ds, dt) = /(p) ds dt, p = (s^-|-t^)^/^. Assume further that P^f{p) dp < 

00 , for 1 < A: < 5. Our aim is to give an explicit upper bound for the density 
of the probability distribution of M/ where I is the unit disc. Using (9) 
which is a consequence of Theorem 3.1 and the invariance of the law of the 
process, we have 

Fp^x) <7rE(P(X"(0,0))/A(0,0) =«,X'(0,0) = (0,0)) 

X PX(0,0),A'(0,0)('*^) (0)0)) 

( 21 ) 

+ 27rE(P(Apl,0))/X(l,0) =xx,X'(l,0) =0)p^(^,,)_^,(^_„)(xx,0) 

= h + l2. 

We denote by X, X', X" the value of the different processes at some point 
(s,t); by the entries of the matrix X"; and by and 4* the 

standard normal density and distribution. 

One can easily check that: X' is independent of X and X", and has vari¬ 
ance nJ^Id] Xg^ is independent oi X, X' X'J^ and X[[, and has variance ^ J5. 
Conditionally on A = u, the random variables X'J^ and X^[ have 


expectation: 
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We obtain 


variance: 

covariance: 


— J5 - (vrJs)^; 


h 



—J5 - (vrJa)^ 


1/2 


ip{bu) + 'Kj3U^{bu) 


with b = (3^/4jgr^(^j3)2)i/2 ■ As -^1 we remark that, conditionally on X = u, 
X'Jg + Xj.^ and X'J^ — X"^ are independent, so that a direct computation gives 


h = 


( 22 ) 


Svr 


■(/?(u)E 


(crr^l 2ttJ^u) ^ (r /2 + %)l{ar;i<27rJ3M} 


^ ^{(Q:7?l-27rJ3ti)2-^(77|+77|)>0} 


where r]i,rj2, rjs are standard independent normal random variables and = 
2 ttJ^ — J3 . Finally we get 





+ a^ 


(?x^)^{a — cx) 


+ [2aQ: — o?{a — cx)]ip{a — ca:)]x(/?(x) dx, 


with a = 27 rJ^u, c 



4. First derivative, second form. We choose, once for this entire section, 
a finite atlas A for I. Then, to every t £ I it is possible to associate a fixed 
chart that will be denoted When t £ dl, (pt{Ut) can be chosen to 

be a half ball with belonging to the hyperplane limiting this half ball. 
For t £ /, let Vt be an open neighborhood of t whose closure is included in 
Ut and let t/’t be a function such that on 14 ; V’t = 0 on U^. 

1. For every t€ I and s £ I we define the normalization n{t,s) in the fol¬ 
lowing way: 

(a) For s £ Vt, we set “in the chart” s) = ^||s — t|p. By 

“in the chart” we mean that ||s — 1\\ is in fact ||(/>t(t) — ())t('S)||- 

(b) For general s, we set n{t,s) ='(/^t(s)ni(t, s) -|- (1 —'iptis))- 

Note that in the flat case, when the dimension d of the manifold is equal 
to the dimension N of the ambient space, the simpler definition n{t, s) = 
^||s — works. 
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2. For every t e dl and s G /, we set ni{t, s) = |(s — t)Ar| + — t|P, where 

{s — t)]sf is the normal component of (s — t) with respect to the hyperplane 
delimiting the half ball (ptiUt)- The rest of the definition is the same. 

Definition 4.1. We will say that / is helix-function—or an /i-function— 
on I with pole t G / satisfying hypothesis k integer, A: > 1, if: 

(i) / is a bounded function on I \ {t}. 

(ii) f{s) := n{t,s)f{s) can be prolonged as function of class on I. 

Definition 4.2. In the same way X is called an h-process with pole 
t G I satisfying hypothesis k integer, A: > 1, if: 

(i) Z is a Gaussian process with paths on /\ {t}. 

(ii) For t G /, Z_{s) := n{t,s)Z{s) can be prolonged as a process of class 
on /, with Z(t) = 0, Z'(t) = 0, If si,..., Sm are pairwise different points 

of I \ {i}, then the distribution of (t), ■ ■ ■, ^ ^ 

does not degenerate. 

(hi) For t G dl; Z(s) := n{t, s)Z(s) can be prolonged as a process of class 
on I with Z(t) = 0, Z (t) = 0, and if si,..., Sm are pairwise different points 
of / \ {t}, then the distribution of Z)y(A), Z^^^(t),..., Z^^)(t), Z(si),..., Z(^^(si),... ,Z^^\sm) 
does not degenerate. Z^j^{t) is the derivative normal to the boundary of I at 

t. 


We use the terms “h-function” and “/i-process” since the function and 
the paths of the process need not extend to a continuous function at the 
point t. However, the definition implies the existence of radial limits at t. 
So the process may take the form of a helix around t. 

Lemma 4.1. Let X be a process satisfying (Hk,k >2), and let f be a 
function I ^TZ. 

(a) For t G I, set for s G I ,s f^t, 

X{s) = a\X{t) + {bi,X'{t)) + n{t, s)X\s), 

where o* and 6* are the regression coefficients. 

In the same way, set 

f{s) = alf{t) + {blJ'{t)) + n{t, s)f\s), 

using the regression coefficients associated to X. 

(b) For t G dl, s gT,s f^t, set 

X{s) = dixit) + iblX'it)) + nit, s)X\s) 


and 


fis) = difit) -h {hi, fit)) -h nit, s)fis). 
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Then s X^{s) and s /*(s) are, respectively, an h-process and an 
h-function with pole t satisfying 


Proof. We give the proof in the case t € I, the other one being similar. 
In fact, the quantity denoted by 2f^{s) is just X(s) — a\X{t) — {h\,X'{t)). 
On L^(n,P), let n be the projector on the orthogonal complement to the 
subspace generated by X(t),X'{t). Using a Taylor expansion, 


X{s)=X{t) + {{s-t),X'{t)) 

+ II^“'S|P / X" {{1 — a)t + as)[v,v\{l — a) da 
Jo 

with V = ||^~*|| . This implies that 

' [ X"— a)t + as)[v,v]{l — a) da 
Jo 


(23) z(s)=n 




which gives the result due to the nondegeneracy condition. □ 


We state now an extension of Ylvisaker’s (1968) theorem on the regularity 
of the distribution of the maximum of a Gaussian process which we will use 
in the proof of Theorem 4.2 and which might have some interest in itself. 


Theorem 4.1. Let Z :T be a Gaussian separable process on some 
parameter set T and denote by = sup^,=rp Z(t) which is (a random vari¬ 
able) taking values mPU{+oo}. Assume that there exists do > 0, m_ > —oo 
such that 

m{t) =E{Zt) >m-, cr^(t) = Var(Zt) > do 

for every t ^ T. Then the distribution of the random variable is the 
sum of an atom at +oo and a—possibly defective—probability measure on 
TZ which has a locally bounded density. 

Proof. Suppose first that X:T —> P is a Gaussian separable process 
satisfying Var(W) = Ij E(-^t) > 0, for every t G T. A close look at Ylvisaker’s 
(1968) proof shows that the distribution of the supremum has a density 
Pm^ that satisfies 

(24) pj^x {u) < V’(«) = ^ for every ueTZ. 

In exp(-v^/2)dv 

Let now Z satisfy the hypotheses of the theorem. For given a,b €TZ, a < b, 
choose A G TZ'^ so that |a| < A and consider the process 
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Clearly, for every t gT, 

E(x{t)) = % i“-i+-^ > _ i’"-i+1-1 + i’"-i+> o_ 


a{t) 


ctq 


o-Q 


o-Q 


and Var(X(f)) = 1. So that (24) holds for the process X. 

On the other hand, the statement follows from the inclusion: 


{a<M^<6}c|- 

which implies 


m-\+A |m_|+A b — a 

-^J^-+- 


P{a <M^ <b}< 


<7o cro 

A|m_ \+A) /cTQ+{h-a)/ao 
(|m_|+A)/(To 

V — a + |m_| + A 


= f-V’ 

Ja O'q 


'<p{u) du 
dv. 


o-Q 


□ 


Set now /3(t) = 1. The key point is that, due to regression formulae, under 
the condition {X{t) = u,X'{t) = 0} the event 

Au{X,P) ■.= {X{s)<u,ysGl} 

coincides with the event 

Au{X\p^) := {X\s) < p\s)u,ysG I \ {t}}, 

where X^ and /3* are the h-process and the h-function defined in Lemma 

4.1. 


Theorem 4.2 (First derivative, second form). Let X :I ^TZ be a Gaus¬ 
sian process, I a C°° compact manifold contained in TZ^. Assume that X has 
paths of class and for s ^t the triplet {X{s), X{t), X'(t)) inTZxTZx TZ'^ 
has a nondegenerate distribution. Then, the result of Theorem 3.1 is valid, 
the derivative F[{u) given by relation (10) can he written as 

f;{«) = jE[det{X‘"{t)- 


(25) 


^ Px{t),x'{t){u,A)a{dt) 

+ (-1)'^-^ [ E[det(l*"(t)-^*"(t))nl^„(^y^.)] 
J dl 

^Px(t),X'(t)(^^^)^(dt) 


and this expression is continuous as a function of u. 
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The notation X_ (t) should be understood in the sense that we first define 
and then calculate its second derivative along dl. 

Proof of Theorem 4.2. As a first step, assume that the process X 
satisfies the hypotheses of Theorem 3.1, which are stronger that those in the 
present theorem. 

We prove that the first term in (10) can be rewritten as the first term 
in (25). One can proceed in a similar way with the second term, mutatis 
mutandis. For that purpose, use the remark just before the statement of 
Theorem 4.2 and the fact that under the condition {2f(t) = u,X'{t) = 0}, 
X"{t) is equal to 

Replacing in the conditional expectation in (10) and on account of the 
Gaussianity of the process, we get rid of the conditioning and obtain the 
first term in (25). We now study the continuity of tt -^/(^)- variable 
u appears at three locations: 

(i) in the density Px{t),x'{t){u,0), which is clearly continuous, 

(ii) in 

E[ det (Z*" (t) - (t)n) , 

where it occurs twice: in the first factor and in the indicator function. 

Due to the integrability of the supremum of bounded Gaussian processes, 
it is easy to prove that this expression is continuous as a function of the first 

u. 

As for the u in the indicator function, set 

(26) i,:=dei{X^"{t)-^"{t)v) 

and, for > 0, consider the quantity 
which is equal to 

Apply Schwarz’s inequality to the first term in (27): 

E[e.lA4Xh/3‘)\A„_,(AS/3‘)] < [E(d)P{A4X^/3*)\A„_;,(X^/3*)}]l/2_ 

The event Au{X^,f3^) \ Au-h{X^, can be described as 

Vs e I\{t} ■. X^ [s) — (3^ {s)u <0; 3so G I\{t} :X^{so) — /3^{so){u — h) > 0. 

This implies that P^{so) >0 and that — \\P^\\ooh < sup^ eAfi} X^is)-PKs)u< 
0. Now, observe that our improved version of Ylvisaker’s theorem (Theorem 
4.1) applies to the process s X^{s) — f5^{s)u defined on I\{t}. This implies 
that the first term in (27) tends to zero as h I 0. An analogous argument 
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applies to the second term. Finally, the continuity of Fj{u) follows from the 
fact that one can pass to the limit under the integral sign in (25). 

To complete the proof we still have to show that the added hypotheses 
are in fact unnecessary for the validity of the conclusion. Suppose now that 
the process X satisfies only the hypotheses of the theorem and define 

(28) X%t) = Z,{t)+eY{t), 

where for each e > 0, is a real-valued Gaussian process defined on /, 
measurable with respect to the cr-algebra generated by {X{t):t £ /}, pos¬ 
sessing paths and such that almost surely Zir{t), Z'^{t), Z”{t) converge 
uniformly on I to X{t), X'(t), X"(t), respectively, as e | 0. One standard 
form to construct such an approximation process Z^ is to use a parti¬ 
tion of the unity on I and to approximate locally the composition of a chart 
with the function X by means of a convolution with a C°° kernel. 

In (28), Y denotes the restriction to / of a Gaussian centered stationary 
process satisfying the hypotheses of Proposition 3.1, defined on TZ^, and 
independent of X. Clearly X^ satishes condition {Hk) for every k, since 
it has C°° paths and the independence of both terms in (28) ensures that 
X^ inherits from Y the nondegeneracy condition in Definition 3.1. So, if 
Mj = maxtsi X^{t) and Ff{u) = P{M| < tt}, one has 

Ffiu) = (-l)Y^E[det(X-^"(f) 


(29) 


X Px^{t),x^'{t){u,0)a{dt) 

_l_ (_ 1 )'^ 1 /" E[det(X (t)—^ {t)u)lA^{X‘^k 0 ^*)] 

J dl 


We want to pass to the limit as e J, 0 in (29). We prove that the right-hand 
side is bounded if e is small enough and converges to a continuous function 
of tt as e I 0. Since Mf Mj, this implies that the limit is continuous and 
coincides with F[{u) by a standard argument on convergence of densities. 
We consider only the first term in (29); the second is similar. 

The convergence of X^ and its first and second derivative, together with 
the nondegeneracy hypothesis, imply that uniformly on t £ I, as e i 0, 
PA®h),x®'(t)(^)0) The same kind of argument can be used 

for det{ 2 F^''{t) — on account of the form of the regression coeffi¬ 

cients and the definitions of X^ and /3*. The only difficulty is to prove that, 
for fixed u, 
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where = Au{X^\C = Au{X\(5^). 

We prove that 

(31) a.s. ICe ^ Ic as e I 0, 

which implies (30). First of all, note that the event 

L = | sup (X*(s) —/3*(s)tt) = 0 

UeAhl 

has zero probability, as already mentioned. Second, from the definition of 
X^{s) and the hypothesis, it follows that, as e | 0, X‘^'*{s), converge 

to X^{s),(5^{s) uniformly on I \ {t}. Now, if cu ^ C, there exists s = s(a;) G 
I \ {f} such that X^{s) — f3^{s)u > 0 and for e > 0 small enough, one has 
X^\s) — I3^^{s)u > 0, which implies that uj ^Ce- 
On the other hand, let a; G C \ L. This implies that 

sup {X^{s) — f3^{s)u) < 0. 

s£l\{t} 

From the above-mentioned uniform convergence, it follows that if e > 0 is 
small enough, then sup^g/\{i|(W^*(s) — (5^*’{s)u) < 0, hence uj G Cg. Equation 

(31) follows. 

So, we have proved that the limit as e | 0 of the first term in (29) is equal 
to the first term in (25). 

It remains only to prove that the first term in (25) is a continuous 
function of u. For this purpose, it suffices to show that the function u ^ 
V{Au{X^, (5^)} is continuous. This is a consequence of the inequality 

\P{A^+h{X\(5^)} - V{Au{X\l3^)}\ 

<\h\ sup |/3‘(s)| 

and of Theorem 4.1, applied once again to the process s ^ X^(s) — P*{s)u 
defined on I \ {t}. □ 

5. Asymptotic expansion of F'{u) for large u. 

Corollary 5.1. Suppose that the process X satisfies the conditions of 
Theorem 4.2 and that in addition E(At) = 0 and Var(Xi) = 1. 

Then, as -foo, F'{u) is equivalent to 

( 32 ) 

where A(t) is the variance-covariance matrix of X'{t). 


<P 


sup (X*(s) — /3*(s)u) 
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Note that (32) is in fact the derivative of the bound for the distribution 
function that can be obtained by Rice’s method [Aza’is and Delmas (2002)] 
or by the expected Euler characteristic method [Taylor, Takemura and Adler 
(2004)]. 

Proof of Corollary 5.1. Set r(s,t) := E(A(s), A(t)), andfori,j = l,d, 
:= —r{s,t), 

OSi 

Eor every t,i and j, = 0, Aij{t) = Thus X{t) 

and X'{t) are independent. Regression formulae imply that a* = r{s,t), 
P^{s) = This implies that /3*(t) = A(t) and that the possible lim¬ 

its values of /3*(s) as s — > t are in the set {v^K{t)v:v G 5"’*“^}. Due to the 
nondegeneracy condition these quantities are minorized by a positive con¬ 
stant. On the other hand, for s^t, f3^{s) > 0. This shows that for every t € I 
one has inf^gj/3*(s) > 0. Since for every tG I the process A* is bounded, it 
follows that a.s. lyi„(xh/3*) ^ 1 as u ^ -|-oo. Also 

det {2L^''{t) — l3^"{t)u) ~ (—l)‘^det(A(t))u‘^. 

Dominated convergence shows that the first term in (25) is equivalent to 

J^u'^det(A*)(27r)-^/2g-nV2(27r)-'^/2(det(A*))-^/^dt 

The same kind of argument shows that the second term is 0(tt‘’*“^e““^/^), 
which completes the proof. □ 
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